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Abstract. It is hereby established that, in Euclidean spaces of finite dimension, bounded self- 
contracted curves have finite length. This extends the main result of [6: concerning continuous 
planar self-contracted curves to any dimension, and dispenses entirely with the continuity re- 
quirement. The proof borrows heavily from a geometric idea of [13] employed for the study 
of regular enough curves, and can be seen as a nonsmooth adaptation of the latter, albeit a 
nontrivial one. Applications to continuous and discrete dynamical systems are discussed: con- 
tinuous self-contracted curves appear as generalized solutions of nonsmooth convex foliation 
systems, recovering a hidden regularity after reparameterization, as consequence of our main 
result. In the discrete case, proximal sequences (obtained through implicit discretization of a 
gradient system) give rise to polygonal self-contracted curves. This yields a straightforward 
proof for the convergence of the exact proximal algorithm, under any choice of parameters. 
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1. Introduction 

1.1. Motivation and state-of-the-art. Self-contracted curves were introduced in [6l Defini- 
tion 1.2.] to provide a unified framework for the study of convex and quasiconvex gradient 
dynamical systems. Given a possibly unbounded interval I of K, a map 7 : I — > W 1 is called 
self- contracted curve, if for every [a, b] C /, the real-valued function 

t€ [a,fe]^d( 7 (t), 7 (6)) 

is non-increasing. This notion is purely metric and does not require any prior smoothness/continuity 
assumption on 7. 

So far, self-contracted curves are considered in a Euclidean framework. In particular, given a 
smooth function / : IR n -4 M, any solution 7 of the gradient system 

= -V/( 7 (*)) * > 0, 
= x G M n 

is a (smooth) self-contracted curve, provided that / is quasiconvex, that is, its sublevel sets 

[/ < 0\ := {x € R n : f(x) < (3} (f3 G K) 

are convex subsets of M n [51 Proposition 6.2]. Self-contracted curves also appear in subgradient 
systems, defined by a (nonsmooth) convex function /, see [6j Proposition 6.4]. In this case, the 




(1.1) 



first equation in (1.1) becomes the following differential inclusion 

Y(t) € -df(a(t)) a.e., (1.2) 

where a.e. stands for "almost everywhere" and the solutions are absolutely continuous curves 
(see [2] for a general theory). 

A central question of the asymptotic theory of a general gradient dynamical system of the form 



( 1.1 ) is whether or not bounded orbits are of finite length, which if true, yields in particular their 
convergence. This property fails for C°° smooth functions ([EJ p. 12]), but holds for analytic 
gradient systems ([H]), or more generally, for systems defined by an o-minimal (tame) function 
(PI). In these concrete estimation of the length of (sub)gradient curves is obtained in 

terms of a so-called Lojasiewicz type inequality, intrinsically linked to the potential /, see pQ for 
a survey, which also includes extensions of the theory to subgradient systems. Notice however 
that convex functions fail to satisfy such inequality, see [TJ Section 4.3] for a counterexample. 

In [13] a certain class of Lipschitz curves has been introduced (with no specific name) to 
capture the behaviour of orbits of quasiconvex potentials. Unlike self-contracted curves, this 
notion makes sense only in Euclidean spaces (it uses orthogonality) and requires the curve to be 
Lipschitz continuous. The main result of [13] asserts that the length of such curves is bounded by 
the mean width of their convex hull, a fortiori, by the mean width of any convex set containing 
the curve. Recently, the authors of [7] extended the result of [13] to 2-dimensional surfaces 
of constant curvature, naming these curves as (G)-orbits. We shall prefer to call these curves 



self-expanding, see Definition 2.2 below. The choice of this terminology will become clear in 
Section O 



Formally the result, as announced in |13l IX], applies only under the prior requirement that 
the length of the curve is finite, since the proof is given for self-expanded curves parameterized by 
the arc-length parametrization in a compact interval. This restriction is removed in Section [2.1| 



via a simple continuity argument (see Corollary 2.4). As a result, the smooth orbits of (|1.1 



for a quasiconvex potential, as well as the absolutely continuous orbits of (1.2) for a convex one 
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have finite length. In both cases the bound depends only on the diameter of the initial sublevel 
set. 

The work [13] was unknown to the authors of [6], who tackled the same problem in terms 
of the aforementioned notion of self-contracted curve. The definition of self-contracted curve 
does not require any regularity neither on the space nor on the curve. In particular, such 
curves need not be continuous, and differentiability may a priori fail at each point. The main 
result of [6j shows that bounded continuous planar self-contracted curves have finite length ([6, 
Theorem 1.3]). This has been used to deduce that in M 2 , smooth orbits of quasiconvex systems 
(respectively, absolutely continuous orbits of nonsmooth convex systems) have finite length. As 
we saw before, this conclusion essentially derives from [13, IX] for any dimension (see comments 
above). Notice however that the main result of [6j cannot be deduced from [13], but can only 
be compared in retrospect, once rectifiability is established. 

On the other hand, not completely surprisingly, Lipschitz continuous self-contracted curves 
and self-expanded curves turn out to be intimately related and can be obtained one from the 



other by means of an adequate reparameterization, reversing orientations (see Lemma 2.8). Re- 
call however that both rectifiability and Lipschitz continuity of the curve are prior requirements 
for the definition of a self-expanded curve, while they are neither requirements nor obvious 
consequences of the definition of a self-contracted curve. 

1.2. Contributions of this work. In this work we establish rectifiability of any self-contracted 
curve in M. n , by extending the result of [6] to any dimension, and to possibly discontinuous curves. 
This is done by adapting the geometrical idea of [13] to the class of self-contracted curves. This 
nonsmooth adaptation is natural but quite involved. Nonsmooth variations of the mean width of 
the closed convex hull of the curve are again used to control the increase of its length, but no prior 
continuity on the parametrization is required and rectifiability is now part of the conclusions. 
Namely, setting T = 7(1) (the image of the curve in ~R n ) and denoting by £(j) its length, we 



establish the following result (see Theorem 3.3) : 



• Every self-contracted curve is rectifiable and satisfies the relation £(7) < Cdiam(r), 
where C > depends only on the dimension. 

In case of continuous curves, the above result allows to consider a Lipschitz reparameterization 



defined by the arc-length, see details in Section 3.2 This leads to the following conclusion: 



• If 7 is a continuous self-contracted curve and T = is bounded, then T is also the 
image of some (Lipschitz) self-expanded curve. 

In particular, the sets of all possible images of continuous self-contracted curves and of self- 
expanded curves coincide. Still, the set of images of all self-contracted curves is much larger (its 
elements are not connected in general). 

In the last two sections, two new applications of self-contracted curves are considered. In 



Section 4.1 we broaden the framework of dynamical systems to encompass nonsmooth convex 
foliation systems, with merely continuous generalized orbits. Limits of backward secants remedy 
the absence of differentiability, leading to a consistent notion of generalized solution in the sense 



of nonsmooth analysis (Definition 4.4). In Theorem 4.6 we show that these generalized solutions 
are self-contracted curves, thus of finite length; in view of the aforementioned result, they can 
also be obtained as "classical" solutions through an adequate Lipschitz reparameterization. 
On the other hand, C 1 smooth convex foliation orbits enjoy a stronger property, the so-called 



strong self-contractedness, see Definition 4.8 and Corollary 4.11 Concerning this latter class 
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we establish in Section 4.2 the following approximation result, with respect to the Hausdorff 



distance, see Proposition 4.13 



Every C 1 -smooth strongly self-contracted curve is a limit of polygonal self-contracted 
curves. 



Finally, in Section 4.3 we provide an elegant application of the notion of self-contracted curve 
in a different framework, that of discrete systems. In particular we establish the following result 



(Theorem 4.17) : 



• Let / be any convex function, bounded from below. Then the exact proximal algorithm 
gives rise to a self-contracted polygonal curve. 

In view of our main result, we obtain a straightforward proof of the convergence of the 
proximal algorithm. The bound over the length of the polygonal curve yields a sharp estimation 
on the rate of convergence, which appears to be entirely new. Notice that the convergence result 
is independent of the choice of the parameters. 
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7) and Laboratory Jacques Louis Lions. The stay was supported by the program "Research in 
Paris" offered by the Ville de Paris (Mairie de Paris). This author thanks the host institution 
and Ville de Paris for its hospitality. 



2. Notation and preliminaries 

Let (W n , d, Jz? n ) denote the n-dimensional Euclidean space endowed with the Euclidean dis- 
tance d(x,y) = \\x — y\\, the scalar product (•,•), and the Lebesgue measure Jzf n . We denote 
by B{x,r) (respectively, B{x,r)) the open (respectively, closed) ball of radius r > and center 
x £ W 1 . If A is a nonempty subset of W 1 , we denote by §A its cardinality, by conv (A) its convex 
hull and by diam^4 := sup{d(x,y) : x,y £ A} its diameter. We also denote by 'mt(A), A and 
dA the interior, the closure and respectively, the boundary of the set A. 

Let now K be a nonempty closed convex subset of W 1 and no £ K. The normal cone Nk(uo) 
is defined as follows: 

N K (u ) = {d £ 1™ : (v,u-u ) < 0,Vu € K}. (2.1) 

Notice that Nk(uq) is always a closed convex cone. Notice further that «o £ K is the projection 
onto K of all elements of the form uq + tv, where t > and v S Nk(uq). 

The Hausdorff distance between two nonempty closed subsets K i, K 2 of M n is given by the 
formula 

d H (K 1 ,K 2 ) := inf le > : K 2 C (J B(z, e) and K x C |J B(z, e) \. 

^ zeKi zGK 2 ' 

A sequence of closed sets {Kj}j in W 1 is said to converge with respect to the Hausdorff distance 
to a closed set K C W 1 if d H (Kj,K) ->■ as j — > 00. 
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Throughout the manuscript, / will denote a possibly unbounded interval of M. In this work, a 
usual choice for the interval will be I = [0, Tqq) where T^ 6lU {+00}. A mapping 7 : / — > M. n 
is referred in the sequel as a curve. Although the usual definition of a curve comes along with 
continuity and injectivity requirements for the map 7, we do not make these prior assumptions 
here. By the term continuous (respectively, absolutely continuous, Lipschitz, smooth) curve we 
shall refer to the corresponding properties of the mapping 7 : I — > W n . A curve 7 is said to be 
bounded if its image, denoted by T = 7(1), is a bounded set of 1". 

The length of a curve 7 : I — > M. n is defined as 

m— 1 

£(7) :=sup{ J>fr(*i).7(ti+i))}, (2.2) 

i=0 

where the supremum is taken over all finite increasing sequences to < t\ < ■ ■ ■ < t m that lie in 
the interval /. Notice that ^(7) corresponds to the total variation of the function 7 : / — > W 1 . 
Let us mention for completeness that the length £(7) of a continuous injective curve 7 is equal 
to the unidimensional Hausdorff measure l-L l {T) of its image, see e.g. Th.2.6.2], but it is 
in general greater for noncontinuous curves. In particular we emphasis that for a piecewise 
continuous curve, the quantity £(j) is strictly greater than the sum of the lengths of each pieces 
but we still call £(7) the length of 7. A curve is called rectifiable, if it has locally bounded length. 

Let us finally define the width of a (nonempty) convex subset K of W 1 at the direction u G S n_1 
as being the length of its orthogonal projection P U (K) on the 1-dimensional space Mu generated 
by u. The following definition will play a key role in our main result. 

Definition 2.1 (Mean width). The mean width of a nonempty convex set K C M n is given by 
the formula 

W(K) = — f J£\P u {K))du, 
where du denotes the standard volume form on and 

n7r n/2 

~r(f + i)- 

2.1. Self-expanded curves. Let us now recall from [133 the definition and the basic properties 
of a favorable class of Lipschitz curves, which has been studied thereby without a specific name. 
In the sequel we call these curves self-expanded. 

Definition 2.2 (Self-expanded curve). A Lipschitz curve 7 : I — > W 1 is called self-expanded 
curve if for every t £ I such that 7'(i) exists, we have that {j ; (t), 7(t) — 7(1*)) > for all u £ I 
such that u < t. 

In |13| 3. IX.] the following result has been established concerning self-expanded curves. 

Theorem 2.3 ([13, 3. IX.]). Let 7 : 1 — > W 1 be a self-expanded curve of finite length. Then there 
exists a constant C > depending only on the dimension n such that 

<C diamK, (2.3) 

where K is any compact set containing V = 7(1). 

Notice that, formally, the above result requires the curve to have finite length. Nevertheless, 
the following limiting argument allows to obtain a more general conclusion for bounded self- 
expanded curves. 
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Corollary 2.4 (Bounded self-expanded curves have finite length). Every bounded self-expanded 



curve 7 : I — > M. n has finite length and (2.3) holds 



Proof. Since self-expanded curves are rectifiable by definition, and because reparameterizing 
7 does not change the statement, we may assume that 7 is parameterized by its arc- length on 
I = [0,^(7)). Notice though that in principle £(7) might be infinite. Our aim is precisely to 
show that this is not the case. Indeed, let K be a compact set containing and let {L n } n be 



an increasing sequence of real numbers converging to £(7) 6lU {+00}. Applying Theorem 2.3 
for the curve 7 n : [0, L n ] — > W 1 (restriction of 7 to [0, L n ]), we obtain 

L n = £(ln) < C diamX, for all n > 1. 

This shows that £(7) = lim n _> +00 L n is bounded and satisfies the same estimate. □ 

2.2. Self-contracted versus self-expanded curves. The aim of this section is to prove that 
Lipschitz continuous self-contracted curves and self-expanded curves give rise to the same im- 
ages. Moreover, each of these curves can be obtained from the other upon reparameterization 
(inverting orientation). As a byproduct, bounded Lipschitz self-contracted curves in W 1 have 
finite length. 

Let us recall the definition of a self-contracted curve (see [61 Definition 1.2]). 

Definition 2.5 (Self-contracted curve). A curve 7 : / — > W 1 is called self-contracted, if for every 
t\ < t2 < t% in / we have 

d(7(*i),7(t3))>d(7(*2),7(f3)). (2-4) 
In other words, the function 1 1— > d{^{t), 7(^3)) is nonincreasing on In (—00, £3]. 

Remark 2.6. (i) As we already said before, the definition of self-contracted curve can be given 
in any metric space and does not require any regularity of the curve, such as continuity or 



differentiability. Notice moreover that if 7(ii) = 7(^3) in (2.4) above, then ^(t) = 7(ii) for 
t\ < t < t2] thus if 7 is not locally stationary, then it is injective. 

(ii) It has been proved in [6l Proposition 2.2] that if 7 is self-contracted and bounded, and 
I = [0, Too) with T^ £ K + U {+00}, then 7 converges to some point 700 £ 1" as t — > T c 



(Notice that this conclusion follows also from our main result Theorem 3.3 ) Consequently the 
curve 7 can be extended to 7 = [0,Too]. In particular, if 7 is continuous, then denoting by 
r = 7(1) the image of 7, it follows that the set 

r = 7(i) u { 7 oo} = i{i) 

is a compact subset of W 1 . 

Lemma 2.7 (Property of a differentiability point). Let 7 : I — > W 1 be a self- contracted curve 
and let t be a point of differentiability 0/7. Then 

(j'(t),'y(u) — j(t)) > for all u £ I such that u > t. 

Proof. Assume that 7 is differentiable at t G I, and write j(t + s) = j(t) + sj'(t) + o(s), with 
lim s _;.o s _1 o(s) = 0. Let u 6 I be such that u > t, take s such that < s < u — t and apply 



(2.4) for ti = t, t2 = t + s and £3 = u. We deduce that 

[|7(*)-7(«)ll > IK* + *)- 7(011- 
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Since 7 (t + s) — j(u) = 7(f) — 7(14) + sj'(t) + o(s), substituting this in the above inequality and 
squaring yields 

> [| 7 ( t + a )_ 7 („)||2_|| 7 (t)_ 7 ( u )[|2 

= 2(«V (t) + o(s), 7 (f) - 7 («)) + ||aV (a) + o(a) || 2 
= 2 S (7 / (i) )7 (i)-7(u)) + o( S ). 
Dividing by s, and taking the limit as s tends to + we get the desired result. □ 

Given a curve 7 : I 1— > M n , we denote by I~ = —I = {— t ; t £ 1} the opposite interval, and 
define the reverse parametrization j~ : I~ >— y R n of 7 by 7~(f) = 7 (— f) for t € I - . 

Lemma 2.8 (Lipschitz self-contracted versus self-expanded curves). Let 7 : J — > W 1 be a 

Lipschitz curve. Then 7 is self-contracted if and only if '7" is a self-expanded curve. 



Proof. If 7 : / — > M. n is Lipschitz and self-contracted, then Lemma 2.7 applies, yielding directly 
that 7~ is a self-expanded curve. Conversely, suppose that 7" is a self-expanded curve. This 
means that 7 is Lipschitz and (after reversing the orientation) that 

(7' (f), 7(14) - 7(f)) > for u € 7 such that u > t (2.5) 

whenever 7 '(i) exists and is different from 0. 

Fix now any t^ £ I and define the function 

/(*) = \h(t) - 7(t 3 )|| 2 , for alU G /, t < t 3 . 

By Rademacher's theorem the Lipschitz continuous functions 7 is differentiable 1 -almost ev- 
erywhere, and so / is too, and /'(f) = (-y' (t),~/(t) — 7^3)) for almost all t G ID (—00,23). 

If i < t 3 and 7 ; (t) / then (|2j5} above (for « = t 3 ) yields that /'(t) < 0. Otherwise, 
/'(t) = 0. It follows that the Lipschitz function / satisfies /'(f) < Jzf^almost everywhere, 



thus it is nondecreasing. This establishes (2.4). Since t 3 has been chosen arbitrarily, the proof 



is complete. □ 

3. RECTIFIABILITY OF SELF-CONTRACTED CURVES 

In [6l Theorem 1.3] it has been established that bounded self-contracted continuous planar 
curves 7 : [0, +00) — > M 2 have finite length. In this section we improve this result by dropping 
the continuity assumption, and we extend it to any dimension. Precisely, we establish that the 
length of any self-contracted curve 7 : I — > W 1 lying inside a compact set is bounded by a 
quantity depending only on the dimension of the space and the diameter of the compact set, see 



the forthcoming Theorem 3.3 



3.1. Proof of the main result. The proof makes use of the following technical facts. 

Lemma 3.1 (Saturating the sphere). Let S C S™" 1 be such that (x,y) < 1/2 for all x,y G 
E, x 7^ y. Then S is finite and jJS < 3 n . 

Proof. For any x, y G S, x 7^ y, we have \\x — y|| 2 = 2 — 2(x,y) > 1. Therefore, the open balls 
{B(x, l/2)} x£ z are disjoint and they are all contained in the ball B(0, §). Set u n = 5? n {B(0, 1)) 
(the measure of the unit ball) ; then 

(8E) co n = sr{ (J 1/2)) < JST (B(0,3/2)) = o; n Q)" 



A. DANIILIDIS, G. DAVID, E. DURAND-CARTAGENA, A. LEMENANT 



and so ftS < 3". 

Lemma 3.2 (Hemisphere lemma). Let S C S™" 1 be a set satisfying 



(x,y) > 



Then there exists Q S S n suc/i i/iai 



n+l 



2n+l 



/or a// x, y € E. 



for all iGS. 



□ 



(3.1) 



Proof. Let {a^jig/ be a family of points in S, maximal with respect to the property that 

(xi,Xj) < - for all i ^ j in I . (3-2) 

(Notice that such a family can easily be constructed by induction.) Applying Lemma |3.1| we 
deduce that m := JJI < 3™. Set 

t> := y^Xj 
is/ 

and let y 6 S be an arbitrary point. If y does not belong to the family {xj}, g / then by maximality 
of the latter, there exists some io such that (xi ,y) > |. If on the contrary, y = xi for some i, 
then we take io = i. In view of (3.1) we obtain 



(v,y) = (x io ,y) + ^T{xi,y) > - - (m - 1) f-J 



n+l 



(3.3) 



whence, recalling that m < 3 n , we get (y, v) > 3 ( n+1 ). This shows in particular that u / 0. 
Let us now compute \\v\\. We write 

II 1 1 2 

IMI 2 = E 3 * =^2\\xi\\ 2 + E = m + E (3.4) 



where 



Using (3.2) we deduce that | E \ < m(m — l)/2. Therefore (3.4) yields 

m(m — 1) m 



\v\\ <m-\- 



So dividing by \\v\\ in (3.3) and setting 



we obtain 



(C,y) = n — rr , 2/> > 
\\v\\ 



-vm + l)< y (3« + l)<3 2 ™. 



/ 1 / \ n+l 



2n+l 



Since y is arbitrary in E, the proof is complete. 



□ 



We are now ready to prove the main result of this section. 
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Theorem 3.3. Let 7 : / — > W 1 be a self- contracted curve. Then there exists a constant C n 
(depending only on the dimension n) such that 

i(Tf) < C n W(K), (3.5) 

where K is the closed convex hull 0/7(1). In particular, bounded self- contracted curves have 
finite length. 

Proof. The result holds vacuously for unbounded curves (both left-hand and right-hand side 
of (3.5) are equal to +00). Therefore, we focus our attention on bounded self-contracted curves 
and assume that K is compact. We may also clearly assume that n > 2 (the result is trivial in 
the one-dimensional case). 

In the sequel, we denote by T = 7(1) the image of such curve. The set T inherits from / a 
total order as follows: for x,y G T we say that u x is before y " and denote x X y, if there exist 
h,t2 G I, t\ < ti and 7(^1) = x, 7^2) = V- \ix <y and x / y, then the intervals 7 _1 (x) and 
7 _1 (y) do not meet, and for any t\ G 7 _1 (x), ti G 7 _1 (y) we have t\ <ti. In this case we say 
that "x is strictly before y" and we denote x -< y. For x G V we set 

T(x) := {y G T : x X y} 

(the piece of curve after x) and denote by Q(x) the closed convex hull of T(x). 



Claim 1. To establish (3.5) it suffices to find a positive constant e = e(n), depending only on 
the dimension n, such that for any two points x, x' G T with x' X x it holds 

W(J2(a;)) + e||x - < TV(fi(a/)). (3.6) 
Proof of Claim 1. Let us see how we can deduce Theorem |3.3| from the above. To this end, 



let to < t\ . . . < t m be any increasing sequence in /, and set Xi = 7(t«). If (3.6) holds, then 

m—1 m—1 1 m—1 

\h(u + i) - j(k)\\ = 2 \\x i+1 - Xi \\ <-J2 {w(c( Xl ) - w(c(x t+1 ))) 

i=0 i=0 i=0 

= ^(W(C(x )) - W(C(x m ))) < -W{C{x G )) < -W{K), 

since the mean width W(H) is a nondecreasing function of H (the variable H is ordered via the 
set inclusion). Taking the supremum over all choices of to < t± . . . < t m in I we obtain (3.5) for 
C n = l/e. ▲ 



Therefore, the theorem will be proved, if we show that (3.6) holds for some constant e > 
which depends only on the dimension. Before we proceed, we introduce some extra notation. 
Given x,x' inT with x' -< x we set 

xo ■= \x + \x and v := t^- — ^tt. (3.7) 
3 3 1 1 a;' — x\\ 

For the sake of drawing pictures, the reader is invited to think that vq = e\ (the first vector of 
the canonical basis of W 1 ), see Figure 1. Let us also set 

/~ X °|, GS"- 1 , foranyyGr(x). (3.8) 

\\y - x o\\ 

Clearly, xq, vq and Coiu) depend on the points x,x', while the desired constant e does not. 
To determine this constant, we shall again transform the problem into another one (see the 
forthcoming Claim 2). 
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Figure 1. Controlling the tail of a self-contracted curve. 



Claim 2. Let us assume that there exists a constant < 5 < 2 4 , depending only on the 
dimension n, such that for all x, x' in T with x' -< x (and for xq, vq defined by ( |3.7| )), there exists 
v G S"- 1 nB(v , 5) such that 



iMv)) < S 2 for ally GT(x). 



(3.9) 



Then (3.6) holds true (and consequently (3.5) follows). 



Proof of Claim 2. Assume that such a constant 5 and a vector v exist, so that (3.9) holds. 



Set 



V = {v£8 



n-l 



< 



Combining with (3.8) and (3.9) we get 

{v, y — xq) < 0, for all v G V and y G T(x). (3.10) 

Let us first explain intuitively why the above yields (3.6). Indeed, compared with Q(x), fi(x') has 
an extra piece coming from the segment [xo,a;']. This piece is protruding in all directions v G V 
(which are relatively close to vq), while (3.10) bounds uniformly the orthogonal projections of 
Q(x) onto the lines M.v. Therefore, the extra contribution of the segment [xo,x'] in P v (Q(x')) 
becomes perceptible and can be quantified in terms of ||xo — x'\\, uniformly in V. Since the latter 
set V has a positive measure, the estimation (3.6) follows. 

To proceed, observe that Q(x') contains the convex hull of Q(x) U whence 

p v (n( x )) c p v (n(x')). 

In particular, 

U l {P v {Sl(x))) < ^ 1 (P ^ ,(0(x , ))) for v G S n -\ (3.11) 

where P v denotes the orthogonal projection onto the line M.v. Let us now equip the latter with 
the obvious order (stemming from the identification ~Rv = M.) and let us identify P v (xq) with 
the zero element 0. Then (3.10) says that for all directions v in V we have 



Notice that 



supP v (n(x)) < < P v (x') < supP„(fi(x / )). 

... 1 



Fo 



x 



RECTIFIABILITY OF SELF-CONTRACTED CURVES IN THE EUCLIDEAN SPACE AND APPLICATIONS 11 

and that V C B(vq, 5 + 5 2 ). Thus for every v G V we have 

(vq, v) > (vq,vq) — \\vq\\ \\vq — v\\ > 1 — 5 — 5 2 > 7/8. 

This gives a lower bound for the length of the projected segment [xo, x '\ onto M.v, which coincides, 
under the above identification, with P v (x'). Thus 



This yields 



Pv{x) > ^ (-\\X ~ X 



1 



^ ill 

> - \\x — X 
4 



W}{P v {£l(x))) + -\\x - x'\\ < H l (P v (n{x'))) for v € V. 



(3.12) 



Integrating (3.12) for v G V and (3.11) for v G S n 1 \ V, and summing up the resulting 
inequalities we obtain (3.6) with 



e = (4a™)" 1 f 
Jv 



du. 



Notice that this bound only depends on 5, so the claim follows. 



Consequently, our next goal is to determine 5 > so that the assertion of Claim 2 holds. The 



exact value of the parameter S is eventually given in (3.25) and depends only on the dimension 



In particular, it works for any self-contracted curve and any choice of points x,x' 6 T. For the 
remaining part of the proof, it is possible to replace 5 by its precise value. Nevertheless, we 
prefer not to do so, in order to illustrate how this value is obtained. In the sequel, the only prior 
requirement is the bound 5 < 2~ 4 . 



Fix any x, x' in T with x 1 < x and recall the definition of xq, vq in (3.7) and £o(y)> for y G T(x) 
in (3.8). Based on this, we consider the orthogonal decomposition 



iv ®U 



(3.13) 



where U = (vq) is the orthogonal hyperplane to vq. The vector v of the assertion of Claim 2 
will be taken of the form 

v = Ti jcTiP ^3-14) 

where £ is a unit vector in IA. This vector will be determined later on, as an application of 
Lemma 3.2, notice however that for any Q G S n_1 C\IA we get v G S n_1 Pi B(vq,5), as needed. 



We set 



r = {yGr(x); («o,£o(v)><-2a}. 



(3.15) 



Notice that for every y G Tq and any v of the form (3.14) we have 

(«, &(!/)> < "2<5 + IN - u|| < -S < -5 2 . 



Thus (3.9) is satisfied for all y in Tq. 



It remains to choose £ (and adjust the value of 5) so that (3.9) would also hold for y G r(x)\To 



This will be done in five steps. In the sequel we shall make use of the decomposition (3.13) of 
vectors £,o(y), y G T(x), namely: 

£o(l0 = {vo,to(v))vo + tfa(v), 
where ^(y) is the orthogonal projection of £o(y) i n U. 



(3.16) 
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Step 1. We establish that for all y £ T(x) \ {x} it holds 



\x — x 



< 0. 



(3.17) 



6||y - sol 

Indeed, since x' -< x -< y, we deduce from self-contractedness that \\y — x'\\ > \\y — x\\, that is, 
y lies at the same half-space as x defined by the mediatrix hyperplane of the segment [x, x'\. 
Hence, denoting by P VQ the orthogonal projection of W 1 on IR^o (which we brutally identify to 
M. to write inequalities) we observe that 



Pvo(y) < o^o(s) + ^Pv (x) = P V0 (x) + 



1 



x — X 



and consequently 



{v ,y- x Q ) = P VQ (y - xq) < P Vo (x) + -||x - x'\\ - Pv (x ) 



\x — X 



thus, dividing by \\y — xq\\, (3.17) follows 



Step 2. We establish that for all y £ T(x) \ Tq, we have: 

llv-*o||> I ^||s-*'|| 



\(vo^o(y))\ = \My)-^(y)\\<25- 



and 



< ||^ 



< 1. 



(3.18) 
(3.19) 

(3.20) 



Indeed, since y £ T(x) \ To, (3.18) follows easily by combining (3.15) with (3.17). The same 
formulas yield that (vq, £o(y)) G ( — 2(5, 0), whence \{vq, £oO/))| < 25. In view of (3.16), both (3.19) 
and (3.20) follow directly. This ends the proof of Step 2. 



Before we proceed, let us make the following observation, which motivates Step 4 (compare 

(3.21) 



forthcoming inequalities (3.22) and (3.24)). We set 



\y - x\\ 



for all y £ T(x) \ {x}. 



It is easily seen that 

(£(lO,£(*)>>0, for ally, z£T(x)\{x}. (3.22) 
Indeed assuming x -<y -< z (the case y = z is trivial) , the definition of self-contractedness yields 

\\z — x\\ > \\z — y\\ . 

Writing z — y = (z — x) — {y — x) and squaring the above inequality, we obtain 

\\z — x\\ 2 > \\z — x\\ 2 + \\x — y\\ 2 — 2(z — x,y — x) , 



which yields (z — x,y — x) > 0. Dividing by norms gives (3.22) 



Our next objective is to show that vectors £(/(y), y £ T(x) \ Tq satisfy a relaxed inequality of 
type (3.22). We need the following intermediate step. 



Step 3. We show that for every y £ F(x) \ Tq, we have: 

H(y)-Uy)\\ < 325. 



(3.23) 
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Indeed 



[|£(i/)-G> 



< 



y-x 



y-x 



< 



y - x\\ 

\\x — Xq\\ 

\\y - x \\ 

2\\x — xq\ 

\\y - x \\ 



+ 



\y - x\ 
1 



+ 



y-x 



y-xo 



\\y - x\\ 

4 ||x — x' 
3 \\y - x\ 



\y - x \ 
l 

\y - xq\\ 



\\y - x \ 
\y - x \\ 



On the other hand, using (3.18) we get: 

^11 > \\y 



\y 



xo\ 



\%o — x\\ > ( 



1 

125 



3' 



x — x 



165 



Combining the above inequalities, we deduce 

U(y) -&(y)ll < 



165 



1 



thus (3.23) follows thanks to the prior bound 5 < 2 



-4 



4\\x — x 



Step 4- We now establish that for every y, z G T(x) \ Tq, we have: 

(#(y),#(z)>>-65£ 

Let us first calculate the scalar product of the unit vectors £,o(y) and £o(z)- 



(3.24) 



<&(!/), &(*)> = (Z(v)M*)) + My) - £(i/Uo(*)> 

= <£(i/U(*)> + (Z(v)M*)-t(z)) + My) -m,Zo(z)) 
> o - \Mz) - - lieo(y) - C(y)ll > -64 5 , 



where (3.22) and (3.23) have been used for the last estimation. 



Using now the above estimation, the orthogonal decomposition (3.16) recalling (3.19) and the 
prior bound 5 < 2 -4 we obtain successively: 

<#(!/),#(*)> = {to(v)M*)) - <«o,£o(i/)X«o,£o(*)) 

> -645 - 45 2 > -655. 
S'iep 5. We shall now be interested in the set 

s = {^(2/);yer(x)\r }cwn§ n - 1 , 

where ffi(y ) = ffi (y)/[|ffi(y)[| (notice that $(y) ^ from fl3.20P ). In view of ( |3.24| | and the 
estimation (3.20), we obtain that for all y,z £ T(x) \ Tq 

-655 



(%(v),%(z))> 1 _ 
We are now ready to give a precise value for 5, namely, 

■\\3n 



> -815. 



(3.25) 

Recalling that n > 2, under this choice the scalar product between two elements of £ is bounded 
from below by the quantity 

■ 1 \ 3n / 1 \ « 



-815 > -81 



14 



A. DANIILIDIS, G. DAVID, E. DURAND-CARTAGENA, A. LEMENANT 



Therefore, the subset £ of the unit sphere of U ~ W 1 1 satisfies the assumptions of Lemma 3.2 
We deduce that there exists a unit vector (sWfl S n_1 such that 



(C,$(y))>{l) 2n ~\ foraiiyer(*)\r . 



(3.26) 



Taking 5 and ( as above in (3.14), we obtain v G S n n -B(t>o, 5). Let us now verify that (|3.9 



holds. As already mentioned (see comments before Step 1), it remains to consider the case 



y G T(x) \ Tq. Notice that in view of (3.17) 



(v -sc,Uy)) = (vo,Uy))-HC,Uy)) < -HCUy)), 



while in view of (3.13), (3.20) and the orthogonality between vq and £ we get 

In 2n-l 



((My)) = <C,#(i/)) = Uo(y)\\ <C,#(y)) > Vi -&{-)'' 



> 



l\2n 



Assembling the above, and using the fact that \\v - 5(\\ = \[\ + 5 2 > 1 we obtain 

This shows that the assumption made in Claim 2 is always fulfilled. The proof is complete. □ 



3.2. Arc-length reparameterization. Having in hand that the length of a bounded self- 
contracted curve 7 : 1 — > M" is finite (c.f. Theorem 3.3), the curve 7 may be reparameterized by 
its arc-length. This reparameterization is particularly interesting when the curve 7 is continuous. 
In this case, we shall see that the arc-length parametrization is Lipschitz continuous and that 
the image T = 7(1) of the curve can also be realized by a self-expanded curve. 

Set I(t) = I n (—00, t], for t & I, and consider the length function 



L : I 

L(t)- 



>[0,%)] 
*(7|/(t))- 



Thus L(t) is the length of the truncated curve j\i(t) an d consequently the function t h-> L(t) 
is nondecreasing, but not necessarily injective — it is locally constant whenever t 1— > j(t) is 
locally stationary, see also Remark |2.6| (i). Let us assume that 7 is bounded and continuous. 
Then, in view of Remark 



2.6 



(ii), the curve can be extended continuously to I and the set 
T = 7(1) U {700} = 7(1) is a compact connected arc. Moreover, t h-> L(t) is continuous and 
L(7) = [0,^(7)]. It follows easily that L~ 1 (s) is a (possibly trivial) interval [ti , ^2] • Consequently, 
the function 

7 : [0, £(7)] -»• W 1 
7(5) = j(t), for any t G L _1 (s) 

is well-defined. 



(3.27) 



Proposition 3.4 (Hidden regularity of continuous self-contracted curves). Let 7 : I — > W 1 
be a bounded continuous self- contracted curve. Then (3.27) defines a Lipschitz self- contracted 
curve with the same image T = 7(1) = 7([0,-£(7))). In particular, T is also the image of some 
self-expanded curve. 



Proof. Let s, s' G [0, ^(7)] be given. Let t,t' G I be such that s 
that s < s' (thus a fortiori t < t'). We claim that 

L{t')>L(t) + d(j(t), 7 (t')). 



L(t) and s' = L(t'). Assume 



(3.28) 
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Indeed, recalling (2.2), for any e > 0, there exists a finite sequence to < . . . < t m 
satisfying 



in I(t), 



m—l 



Ei =: Y^d(n{ti)MU+x))>L{t)-e. 



i=0 



Since adding an extra point can only make the sum larger, we may assume that t m = t. Then 
add the point t m+ \ = t' and notice that to < . . . < t m < t m+ \ is a sequence in I(t'), thus 



E 

i=0 



d(7(ti),7(ti + l)) 



d( 7 (t), 7(0) <£(*') 



We deduce that 



L(i) - e < Ei = S 2 - d( 7 (t), 7 (0) < MO - d( 7 (t), 7 (0)- 



Taking the limit as e — > and get (3.28). It follows that 

IItV) 



L(t) 



7(«)ll = 117(0- 7(011 <^(0 
which proves that 7 is 1-Lipschitz on [0,^(7)]. 

It follows easily that 7 is self-contracted, while the fact that 7 ([0,^( 7 ))) = 7(1) is straightfor- 



ward. The last assertion follows from Lemma u. 



□ 



Remark 3.5 (Improving the length bound). Although the constants C > in (2.3) (for self- 



expanded curves) and in ( |3.5| ) (for bounded self-contracted curves) stem from the same method, 
the latter is based on nonsmooth arguments encompassing nonregular curves and therefore it 
is not optimized. As a consequence of Proposition 3.4 given a continuous self-contracted curve 
we can use (2.3) to improve a posteriori the constant that bounds the length given in (3.5). 
Notwithstanding, this trick cannot be applied to any bounded self-contracted curve: in the 
discontinuous case, one can still consider the length mapping L(t), and use it to construct a 
larger curve (whose image contains T) with a Lipschitz parametrization, but which is possibly 
not self-contracted. 



4. Applications 

4.1. Orbits of convex foliations. As mentioned in the introduction, self-contracted curves 
are naturally defined in a metric space without prior regularity assumptions. Nonetheless, this 
notion has been conceived to capture the behaviour of orbits of gradient dynamical systems 



for (smooth) quasiconvex potentials in Euclidean spaces (1.1), or more generally, orbits of the 
subgradient semi-flow defined by a nonsmooth convex function, see (1.2). In both cases, orbits 
are at least, absolutely continuous. In this section we broaden the above framework by introduc- 
ing a new concept of generalized solution to a nonsmooth convex foliation, admitting (merely) 
continuous curves as generalized orbits. 

4.1.1. Nonsmooth convex foliations. Let us first give the definition of a convex foliation. 

Definition 4.1 (Convex foliation). A collection {C r } r >o of nonempty convex compact subsets 
of M n is a (global) convex foliation of W 1 if 

n < r 2 C ri C int C r2 (4.1) 

and 

aa = K n \intC . (4.2) 



U 



r>0 
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In view of (4.1), relation (4.2) is equivalent to the fact that for every x 6 1" \ intCb, there 
exists a unique r £ [0, oo) such that x 6 dC r . The above definition is thus equivalent to [6j 
Definition (6.5)]. 

Example 4.2 (Foliation given by a function), (i) The sublevel sets of a proper (coercive) convex 
function / : 1" — > R provide a typical example of a convex foliation. Indeed, if m = min /, we 



set Co = / 1 (m) and C r = [f < m + r\. Then (4.2) follows from the fact that / cannot be 
locally constant outside of Cq. 



(ii) The sublevel sets of a (coercive) quasiconvex function might fail to satisfy (4.2), unless 
int [/ < A] = [/ < A] for all A > min /. This condition is automatically satisfied if the quasiconvex 
function is smooth and every critical point is a global minimizer. A more general condition is 
to assume that / is semistrictly quasiconvex, see for instance [5] for the relevant definition and 
characterizations. 

4.1.2. Generalized solutions and self-contractedness. We consider a possibly unbounded interval 
I = [0, Too) of M, where £ M + U {+00} and a continuous injective curve 7 : I — > W l . For 
every t £ I we define the set of all possible limits of backward secants at 7(1") as follows: 

sec-(r) := !.q£S n ' 1 : q = lim /frj - ^ 1 , (4.3) 

where the notation {t^^ /* r _ indicates that —> r and t^ < t for all A;. The set 

sec + (r) is defined analogously using decreasing sequences {ifcjfc \ T+ ■ The compactness of 
gn-i g uaran tees that both sec~(r) and sec + (r) are nonempty. The following lemma relates 
sec~(r) and sec + (r) with the derivative of 7, in case the latter exists and does not vanish. 

Lemma 4.3 (Secants versus derivative). Assume 7 : I — > W 1 is differentiable at t £ I and 
V(r) / 0. Then 

sec-(r) = (- ,| 7 i r j„ j and sec+(r) - ^ 7 ' (r) 



7'(t)||J v ; U|7'(t) 

Proof. Let q £ sec _ (r) and let {tt}k / be a sequence realizing this limit. Then writing 

l(tk)-l(j) (l(t k ) -7(r)\ / -|r-i fe | 



I|7(<*)-7(t)II V tk-r J\h(tk)-l{r) 

and passing to the limit as {t k }k /* T ~ we obtain the result. The second assertion follows 
similarly. □ 

Assuming j(t) ^ int Cq we denote by r(t) the unique positive number such that 

l(t)£dC r{t) . (4.4) 
We are now ready to give the definition of generalized orbit to a convex foliation. 
Definition 4.4 (Convex foliation orbits). A continuous curve 7 : I — > W 1 is called generalized 



solution (orbit) of the convex foliation {C r } r >o if 7(-0 H int Co = and under the notation (4.3) 
and (|4.4l) we have: 



(i) the function 1 1— > r{t) is decreasing 

(ii) sec-(t) C N Cr{t) (l{t)) for all t £ I. 
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It follows directly from (i) that 7 is injective, and then (ii) makes sense. Also, every orbit of 
a convex foliation is an injective mapping with bounded image (r = 7(1) is contained in the 
compact set C r ( )). Recall the definition of the normal cone Nc r{t) in (2.1 ). If 5C r (t) * s a smooth 
manifold (or more generally, if dimNc, t J r y(t)) = 1) we have sec~(t) = {^ r (t)} where ^ r (t) is the 
unit normal of dC r u\ at j(t). 

In the sequel, we shall need the following lemma, which proof is an easy exercise and will be 
omitted. 

Lemma 4.5 (Criterium for decrease). Let ip : [0, T] — > M be a continuous function, where T > 0. 
Assume that for every < r < T, there exists 5 > such that for all t G (r — 5, r) n [0, T] we 
have 

<p(t) > <p(r). (4.5) 

Then <p is decreasing. 

We now prove that the generalized solutions of a convex foliation are self-contracted curves. 
Notice that in view of Theorem 13.31 this entails that these curves are rectifiable and have a finite 
length. 

Theorem 4.6 (Self-contractedness of generalized orbits). Every convex foliation orbit is a (con- 
tinuous) self- contracted curve. 

Proof. Let 7 : I — > 1" be a continuous injective curve which is a generalized solution of the 



convex foliation {C r }r>o in the sense of Definition 4.4 Since 7(1) C C r iQ\, the curve is obviously 
bounded. Fix t £ I and e > such that r + e G 1. We set 

T e :={7(t) :t >r + e, t G /}. 

Notice that the compact set T e is contained in int C r ( T+e n for all < e' < e. 

Claim 1. There exists a > 0, such that for all v £ -^c r(T ) (7( r )) an d all a; G T £ we have 

|| X ~ 7 S T ||| ) <~a<0. (4.6) 
||x-7(r)|| 



Proof of Claim 1. Assume, towards a contradiction, that (4.6) fails. Since the sets T e and 
Nc r{T) {l{ T )) nS n_1 are compact, we easily deduce that for some unit normal v G Nc^^ir)) 
and 1 £ T e we have (v, x — j(t)) > 0. Since x G int C r t T \, we deduce that for some y G C r i T \ we 
have (v,y — 7(r)) > 0. This contradicts the fact that v G iVc - . (7(7")). Thus (4.6) holds. A 

Claim 2. For every e > there exists J > such that for all t~ ,t + £ I such that t _ G [r — <5, r) 
and t + > r + e we have 

( 7 (t-)_ 7 ( r ), 7 (t+)_ 7 ( r )) < . (4.7) 

Proof of Claim 2. Fix e > and let a > be given by Claim 1. Let N 1 = iV Cr . (T) (7(r)) nS^ 1 
denote the (compact) set of unit normals at 7(1"), let C/ a := A rl +B(0, a) be its a-enlargement, and 
let N a denote the closed convex cone generated by U a . It follows from (4.6 ) that if a > is chosen 
small enough (depending on a), then for all w G A^, and all x G T t , we have (w,x — 7(7*)) < 0. 
Since sec _ (r) C A 1 C U a , we deduce that for some (5 > we have 

7(t) -7(7-) 



[|7(*)-7(r)[| 



G U a , for all t G (r - <5, r) n /. 



This establishes (4.7). 
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Notice that ( gjj ) yields 

h(n-j(t+)\\ > || 7 (r)- 7 (f + )||. (4.8) 

We now prove that the curve 7 is self-contracted. Indeed, fix < t\ < i 2 < £3 in / and consider 
the real- valued function 

f <p:[Q, t 2 ] -> R 

\ p(t) = ||7(t)-7(*3)|| 

Applying Claim 2 for any < r < t 2 and for e = £3 — i 2 we deduce that for some <5i > and 



all t G (t — 5i,t) n [0, t 2 ] we have (/?(i) > </?(r). The conclusion follows from Lemma 4.5 Since 



fi,i 2 ,t3 are arbitrarily chosen we deduce that 7 is self-contracted. □ 



As a consequence of Theorem 4.6 and (BJ Proposition 2.2] the solution curve 7 has finite 



length. In particular, the curve converges as t — > to some limit point 

7oo := lim 7(f) G dC r(Too) 

so that the mapping 7 can be continuously extended to [0, Too], by setting 7(T 00 ) := 7 
Remark 4.7 (Generalized versus classical solutions). Let 7 : I — > W 1 be a convex foliation 



00 • 



orbit and T = Although 7 is merely continuous, Theorem 4.6 guarantees its rectifiability. 

Thus Proposition |3.4| applies, and the curve can be reparameterized by its arc-length to obtain 
a Lipschitz curve with the same image. The new curve 7* : [0,^(7)] —> M. n satisfies ||7*(s)|| = 1 
and 

-7t0) e N c{r(s)) (^(s)) 

for almost all s G [0, 7 (-f)]. It follows that every generalized orbit of a convex foliation gives rise 
to a classical solution (Lipschitz orbit). 

4.1.3. Smooth solutions and strong self-contractedness. Let 7 : J — > W 1 be a continuous curve, 
where I = [0, T^), let T £ I and denote by £l(T) the convex hull of the tail of the curve, that is, 

Q,(T) := conv{7(£) : t > T}. 

Let us give the following definitions. 

Definition 4.8 (Strong self-contractedness). A continuous injective curve 7 : I — > W 1 is called 

• strictly self-contracted, if for ti,tz, £3 in I such that t\ < t 2 < ^3 we have 

Il7(*i)- 7(t3)|| > ||7(*2)-7(*3)||- 

• strongly self-contracted, if for all r G / we have 

sec-(r) CintiVn (r) (7(r)). (4.9) 
Of course every strictly self-contracted curve is self-contracted. A careful look at the proof of 



Theorem 4.6 actually reveals that the orbits of the convex foliation are strictly self-contracted 



curves. The following proposition relates the notions of strict and strong self-contractedness. 

Proposition 4.9 (Strong versus strict self-contractedness). Every strongly self- contracted curve 
is strictly self- contracted. 



RECTIFIABILITY OF SELF-CONTRACTED CURVES IN THE EUCLIDEAN SPACE AND APPLICATIONS 19 



Proof. Fix t\ < t<i < t% in / and apply (4.9) for r = t<i- From the definition of sec (r), see 



(4.3), we deduce that for some 5 > and all t G fa — <£, £2) we have 

(7(*) - 7(<a), 7(*3) - 1(h)) < 0. 



Setting <£>(i) = ||7(i) — 7(^3) ||, the above inequality guarantees that Lemma 4.5 applies. The 



proof is complete. □ 

Our next objective is to show that smooth orbits of a convex foliation are actually strongly 
self-contracted curves. We shall need the following lemma. 

Lemma 4.10 (Differentiability point of a convex foliation orbit). Let 7 : I — > W 1 be a continuous 
orbit of a convex foliation and assume that 7 is differentiate at some r > with 7'(t) 7^ 0. 
Then 

-7 / (r)eintiV n(T) (7(T)). (4.10) 
Proof. In order to simplify notation we assume 7(7") = 0. (There is no loss of generality in 



doing this.) In view of Lemma 4.3 there exist e > 0, such that for all t + G (t, r + e) we have 

{ J^L^tL }> l (4ii) 

Since the compact set 0(r + e) is contained in intC r ( T ), arguing as in Claim 1 of the proof of 



Theorem 4.6, we deduce that for some a > and for alH + > r + e we have 



/_ JZM- ^ ) < - a (4 12) 



Changing a into min(a, 1/2) if necessary, we deduce that (4.12) holds true for all t + > r. We 
deduce that the tail of the curve for t > r is contained in the convex cone 

\x€R n : (/^%,x) <-a\\x\\\, 
I \\r{r)\\ ) 

hence so does Cl(r). This easily yields —"/{t) G int A r Q( T )(7(r)). □ 

Corollary 4.11 (Strong self-contractedness of smooth orbits). Every C 1 convex foliation orbit 
with no stationary point is a strongly self- contracted curve. 

Proof. Parameterizing the curve by its arc-length parametrization we obtain a C 1 curve with 
nonvanishing derivative and the same image V = This new curve satisfies ||7(s)|| = 1, for 



every s G [0,^(7)] and it is also a convex foliation orbit, see Remark 4.7 The result follows by 



combining Lemma |4.3| with Lemma 4.10 □ 



4.2. Polygonal approximations of smooth strongly self-contracted curves. In this sec- 
tion we prove that every strongly self-contracted C 1 curve is a limit of self-contracted polygonal 
curves (with respect to the Hausdorff distance). Let us recall the relevant definition. 



Definition 4.12 (Polygonal approximation). Let 7 : i" — > W 1 be a continuous curve. A polygonal 

m 

line P = |J [-Zfc+i, Zk] is called polygonal approximation of accuracy 5 > for the curve 7, if 
fc=o 

{zk}kC~f(I), z k+ i<z k and d H {P,l{I)) < S. 
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In the above we used the same notation x <y which corresponds to the order on the curve as 



in the beginning of the proof of Theorem 3.3 In the sequel, bounded continuous self-contracted 
curves 7 : / — > W 1 are considered to be extended to /, so that T = 7(1) is a compact set. 

Proposition 4.13 (Self-contracted polygonal approximations). Let 7 : [0,L] M. n be a C l 
strongly self- contracted curve. Then for every 5 > the curve 7 admits a self- contracted poly- 
gonal approximation of accuracy 5 > 0. 



Proof. Since the statement is independent of the parametrization, there is no loss of generality 



in assuming that the curve is parameterized by arc- length. In particular, [I7 (t)\\ = 1 for all 
t G [0, L] where L is the total length of the curve. Moreover, (4.10) holds. 

We shall show that for every 5 > and every To G (0, L] there exists a polygonal approximation 

m 

Ps(T ) 



U [z]c+i, Zk] of the curve j([Tq,L]) of accuracy 5 > 0. Indeed, let To G (0, L] and set 
fc=o 

7(Tq). In view of (4.10) and Lemma 4.3 there exists t~ G [Tq — 5, To) such that 



7(f) - 7 (?o) 
|7(t-)-7(To)|| 



GintiV n(To) (7(To)). 



(4.13) 



Let T\ > be the infimum of all t G [To — 5, To) such that (4.13) holds, and set z\ = 7(2~i). 
Notice that {z\ — zq) G 7Vq(j> )(7(To)), which yields that for any x G fl(To), the function 

1 1->- d(x, z + t(zi - z )) 

is increasing. Thus the curve [21,20] U {^(t) : t > To} is self-contracted. Since t h-> j(t) is the 
length parametrization, it follows that \\zi — zq\\ < |Ti — Tq\ < 5 whence 

Notice finally that 

[zi,2d]Ufi(T )cn(Ti). 

Repeating the above procedure, we build after m iterations, a decreasing sequence To > T\ > 

m 

... > T m+ \ and points Zk = 7(Tfc), such that |J [zk+i, Zk] is a self-contracted polygonal 

fc=0 

approximation of the curve {^(t) : t G [T m+ i,To]} of precision 5 > 0. Notice also that 

m 

O(T )U U [Zk+l,z k ] c n(T m+ i). 

k=0 

The proof will be complete if we show that with this procedure we reach the initial point t = 
in a finite number of iterations, i.e. T m+ \ = for some m > 1. Let us assume this is not the case. 
Then the aforementioned procedure gives a decreasing sequence {T m } m \ T, where T > 0. Set 
z = 7(T) and z m = j(T m ) and notice that {z m } — > z. We may assume that {z m } m > mo C B(z, 5). 
It follows that for every m > mo there exists s m > and x m '■= j(Tm + s m ) such that 



(l 



> 0. 



(4.14) 



(If (4.14) were not true, then we could have taken z m+ \ = z and T m+ i = T a contradiction to 
the definition of T m+ \ as an infimum.) Let us now assume that a subsequence of {x m } m remains 
away of z. Then taking a converging sub-subsequence {^A; m } m and passing to the limit in (4.14) 
as 771 — > 00, we obtain that (j'(T),u) > 0, for some unit vector u in the cone over O(T) — z, 
which contradicts (4.9). It follows that s m — > and {x m } m — > z. 
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For the rest of the proof, let us assume for simplicity that z = and '/(T) = e n = (0, . . . , 1). 
(There is no loss of generality in doing so.) Since 7 is a C 1 curve, it follows that for any e > 0, 
there exists r > such that for all t G (T, T + r) we have 

1 > (e„,7'W) > 1-e (4.15) 

and 

7(0 



Pick any m sufficiently large so that Tfn+Sfh < t, set 2 
Since x,z satisfy (4.14) and (4.16) we deduce that 

< and 



< e. 



(4.16) 

j(Tfn) and x = x m = 7( r m + Sm)- 



(1 



x — z 



hence (e n , 



\z\\ \\x — z\\ 
< e, and the last coordinates satisfy 



x 



z l < e x - z\ 



e ||7(T m + Sfh) - l{Tm)\\ < es t 



On the other hand, using (|4.15|) we deduce 

(e n , gamma (Tm + s))ds > (1 — e)sr 



(4.17) 
(4.18) 



Relations (4.17), (4.18) are incompatible for e < 1/2, leading to a contradiction. 

This proves that for some m we get T m+ \ = and z m+ i = 7(0). The proof is complete. □ 



In view of Corollary 4.11 we thus obtain. 

Corollary 4.14. Every C 1 convex foliation orbit admits self-contracted polygonal approxima- 
tions of arbitrary accuracy. 

4.3. Convergence of the proximal algorithm. The aim of this section is to provide a diffe- 
rent type of application of the notion of self-contracted curves. We recall that for a (nonsmooth) 
convex function / : M n — > R which is bounded from below, an initial point xq G 1" and a 
sequence {U}i C (0, 1], the algorithm defines a sequence {xj}j>o C W 1 called proximal sequence 
according to the iteration scheme, see |16| Definition 1.22] for example: 

Given Xi, define as the (unique) solution of the following (strongly convex, coercive) 

minimization problem 



min I f{x) + — ||a; - Xi\ 

A necessary and sufficient optimality condition for the above problem is 

G df(x i+ i) + tr l (x i+ i - xi) 

or equivalently, 



(4.19) 



•Ei+l Xj 
ti 



G -df(x i+1 ), 



(4.20) 



where df is the Fenchel subdifferential of /. Notice that (4.20) can be seen as an implicit 
discretization of the subgradient system (1.2) or of the gradient system (1.1), in case / is 
smooth. 

The proximal algorithm has been introduced in p3]. We refer to [10] for a geometrical 
interpretation and to [1] for extensions to nonconvex problems. Let us now recall from [TJ 
Section 2.4] the following important facts (we give the short proof for completeness): 
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Lemma 4.15 (Geometrical interpretation). Let {xj}j>o be a proximal sequence for a convex 
function f and set ri := f(x{). Then: 

(i) f(x i+ i) < f(xi) if x i+ i / Xi ; 

(ii) the point Xj+i is the shortest distance projection of Xi on the sublevel set [f < rj+i]. 



Proof. Assertion (i) is obvious, since a?i+i minimizes (4.19), thus for x = Xj we get 

f(x i+ i) + ~ x if < f( x i)- 



Assertion (ii) follows again from (4.19), by considering x G [/ < r^+x]- Indeed: 



/(x i+ i) + — \\x i+ i - Xi\\ 2 < f{x) + — 1| a; - Xi\\ 2 < r i+x + —\\x - Xi\\ 2 . 

Since rj + i = /(xj+i) the above yields 

— Xi\\ < \\x — Xi\\, for all x G [/ < ri+i]. 
The proof is complete. □ 

We obtain the following interesting consequence. 

Proposition 4.16 (Self-contractedness of the proximal algorithm). Let {xi}i>o be a proximal 
sequence for a convex function f. Then the polygonal curve 

P=\J[xi,x i+1 ] (4.21) 

is a self- contracted curve. 

Proof. Let i\ < 12 < 13- Notice first that for every x £ [f < r m+ i], the function 

t l— ^ H^m+l + t\X m — X m -i-i) — x\\ 

is increasing because the sublevel sets are convex. Now Lemma 4. 15| yields that for all i\ < m < 
h, Xi 3 G [/ < f{x m+ i)} and 

ll^m+l x i3 II — ll^m Xi^ ||, 

which inductively yields 

11*^2 II — 11^*1 x iz II* (4.22) 
A simple argument now shows that the polygonal curve is self-contracted. □ 



Combining the above with the length bound given in (2.3) (in view of Remark 3.5) we obtain 
the following convergence result. 

Theorem 4.17 (Convergence of the proximal algorithm). Let f : W 1 — > R be convex and 
bounded from below. Let xq G M n be an initial point for the proximal algorithm with parameters 
{ti}i C (0, 1]. Then the proximal sequence {xi}i converges to some point Xoo G W 1 with a fast 
rate of convergence. That is, there exists c > 0, depending only on the dimension, such that 

^ ^ \\xi X^-flU c||xo -£oo||* 
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Proof. By Proposition 4.16 the polygonal curve P denned in (4.21) is self-contracted. Conse- 



quently, by Remark 3.5 we know that ( |2.3[ ) applies, yielding 



net 



x i+i\\ 5; C n diamK, 

where K is any compact convex set containing In addition, letting i\ = and i% — > +oo 

— II^O -^oo || > 



in (14. 221) we get 



H^i ^oo I 

for all i 6 N, which implies that {x{}i C B(xoo,R) with i? 
for c = 2C„ . 



I^o — %oo\\ and the result follows 

□ 



Theorem 4.17 guarantees the fast convergence of the proximal sequence {xt}i towards a limit 
point, with a convergence rate independent of the choice of the parameters Notice however 
that for an arbitrary choice of parameters, the limit point might not be a critical point of / as 
shows the following example: 

Example 4.18 (Convergence to a noncritical point). Consider the (coercive, C 1 ) convex func- 
tion / : R -> E defined by 



/(*) 



x J , |x| < 1/2 
W-l, |a?| > 1/2 

Then the proximal algorithm {xt}i, initialized at the point xq = 2 and corresponding to the 
parameters t{ = 1/2 J+1 , i E N, converges to the noncritical point x = 1. 



In the above example, the choice of the parameters {ti} in the proximal algorithm has the 
drawback that < +oo. In practice the choice of the step-size parameters {ti}i is crucial to 

obtain the convergence of the sequence {f(xi)}i towards a critical value; a standard choice for 
this is any sequence satisfying ^ tj = +oo, see for instance [5]. 
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